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Abstract 

Group theoretic method for the systematic study of five-quark states with meson-baryon (qq — q 3 ) 
configuration is developed. The calculation of matrix elements of many body Hamiltonian is 
simplified by transforming the physical bases (meson-baryon quark cluster bases) to symmetry 
bases (group chain classified bases), where the fractional parentage expansion method can be used. 
Three quark models, the naive Glashow-Isgur model, Salamanca chiral quark model and quark 
derealization color screening model, are used to show the general applicability of the method and 
general results of constituent quark models for five-quark states are given. The method can also 
be useful in the calculation of meson-baryon scattering and the study of the five-quark components 
effect in baryon structure. The physical contents of different model configurations for the same 
multi-quark system can also be compared through the transformation between different physical 
bases to the same set of symmetry bases. 

PACS numbers: 12.39.Mk, 12.39.Jh 
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I. INTRODUCTION 



After 40 years of quark model study, the idea about baryon and meson is about to 
go beyond the naive picture: baryon g 3 and meson qq. The proton spin puzzle could be 
explained by introducing the q 3 qq component in the quark model |l| . Recently to understand 
ihe strange magnetic moment of proton, the five-quark component of proton was proposed 
2]. The recent progresses on meson spectroscopy called for the tetra-quark components 
3|. Although the pentaquark states claimed by experimental groups few years ago might 
be disappeared and the multi-quark states might be hard to be identified, the multi-quark 
study is indispensable for understanding the low energy quantum chromodynamics (QCD), 
because the multi-quark states can provide information unavailable for qq meson and q 3 
baryon, especially the property of hidden color structure. Generally a multi-body interaction 
multi-channel coupling calculation is needed in the multi-quark study. Therefore a powerful 
method is necessary. The fractional parentage (fp) expansion technique has been proven as 
a powerful method for few-body problem. 

To full play the powerful group theoretic method for quark model calculation of multi- 
quark system, one needs not only the fractional parentage expansion coefficients of these 
multi-quark states but also a relation between various quark model states (hereafter called 
physical bases) and the group theoretical classification states (hereafter called symmetry 
bases). Such a method has been developed and successfully applied in systematic search 
of dibaryon [4] and five-quark states with Jaffe-Wilczek di-quark configuration Q], where 
the physical bases were transformed to symmetry bases first, then the many-body matrix 
elements calculation of Hamiltonian (with two body interaction) on the symmetry bases was 
done by means of fractional parentage expansion, i.e., the many-body matrix elements can 
be reduced to overlap and two-body matrix elements. At last the matrix elements on the 
symmetry bases were transformed back to physical bases. 

The present work is to provide the transformation coefficients between physical bases and 
symmetry bases of five-quark systems to facilitate the calculation of many body Hamiltonian 
matrix elements. The physical bases discussed in this paper is the the meson-baryon bases. 
For Jaffe-Wilczek (JW) di-quark bases, see the first paper of this series ({5|. The method is 
applied to the quark model calculation of five-quark systems with three quark models: the 
naive quark model, i.e., the Glashow-Isgur model [6|, the Salamanca version of chiral quark 
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model [3] and the quark derealization color screening model (QDCSM) {4,8]. The general 
applicability of this method is verified through these model calculations. 

Even though the multi-quark system might be hard to be identified, the method is still 
useful in the study of hadron-hadron scattering and the multi-quark components effect 
in hadron in the framework of quark model. The transformation between physical bases 
and symmetry bases is also useful in the study of the physical contents of different model 
approaches for the same multi-quark system. This report gave an example where the different 
model approaches in the study of pentaquark states are compared. 

In section II, the physical bases and symmetry bases are introduced and the transforma- 
tion between them is derived. The fractional parentage technique applied to calculate the 
matrix elements on the symmetry bases is also explained in this section. Section III explains 
three quark models we used. The results of the systematic calculation of five-quark system 
in the u, d, s 3-flavor world are given in section IV. The last section gives the summary. 

II. PHYSICAL BASES AND SYMMETRY BASES 

The physical bases are constructed as follows, first the wave function of each quark cluster 
is constructed based on the group chain classification 

[1»] M [u\ [c] [fj] [f] I Y J 
SU 36 D SU* x {sUig D SU^ x [SU 6 D (SU{ D SU^ x Uf) x SU^] } (1) 

(the Young diagrams or quantum numbers for each group are also shown in Eq.flTJ), then 
the quark cluster wave functions of the system was coupled to definite color, spin and isospin 
quantum numbers by Clebsch-Gordan (CG) coefficients of color SU3, spin SUg and isospin 
SU2 group, and finally anti-symmetrized. For meson-baryon {qq — q 3 ) configuration, the five 
quarks are separated into two clusters: the meson and baryon clusters with separation S. 
The physical basis (meson-baryon basis) is defined as 

^ak(q A q) = A [^3(qiq2q3)i>2(q4q~5)\w I c M I Mj > ( 2 ) 

where A is a normalized antisymmetrization operator for four quarks. ^3 and ip2 are wave- 
functions of baryon and meson, respectively. [ ] means coupling in terms of the SU3, SU2, 
SU2 CG coefficients so that it has color symmetry [c]W c , isospin IMj, and spin JMj. 



a = (YIJ), k represents the quantum numbers z/ 3 , v 2 , C3, • • • , J2, which specify the baryon 
and meson states. In order to relate the physical bases to the symmetry bases, the 4-quark 
cluster basis is introduced 



4 M l4 Mj 4 > 



(3) 



where a 4 = (Y4/4J4), k 4 represents the quantum numbers is 3 , v\, c 3 , • • • , J\. Coupling the 
anti-quark basis to 4-quark basis Eq.(j3]) in terms of SU3, SU2, SUg CG coefficients, gives 
the 5-quark basis 



[c]IJ 
WMjMj 



It relates to the physical basis through the SU3, SU2, SUg Racah coefficients, 

^ak(q 4 q) = U {czCxcci; c^U {hIiIIx,IJ 2 )U{JzJiJJx, J4J2) 
x *a 4 fc 4 (? ))WWJVk 

where U's are Racah coefficients, defined as, 



WMiMj 



(4) 



C/(c 3 CiCCi;c4C 2 ) 



[c 4 ]Wc 



c]W c 



c, 



c 2 ]W c . 



a 



[c]W c 



all w 



The Racah coefficients can be found in Ref. 9J. 

The symmetry basis of 4-quark system can be defined as, 



[c^WMiMYJMJ^Mj, 



(5) 



which is the basis vector belonging to the irreducible representations of group chain Eq.([T]) 
with n = 4. In Eq.Q K4 stands for z/ 4 /i 4 / 4 . Coupling the antiquark to Eq.flSJ), the symmetry 
bases for 5-quark is obtained 

[c]IJ 



®aK{q q) = ®a 4 K 4 (q ))^Ij(gs) 



WM!Mj 



(6) 



To relate the symmetry basis to physical basis, first we express the four-quark cluster 
basis Eq.flSj) in terms of the symmetry basis Eq.([5]) 1(3, [ll|, \\\ . 



u 4 [i 4 f 4 



,[h]Y 4 I 4 
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C's are the isoscalar factors (ISFs) of SU mn DSU m xSU n , which can be obtained from Chen's 
book [9j. Then, the physical bases and symmetry bases for 5-quark can be transformed to 
each other by 

K 

= {UUU} Racah {CCC} lSF ^ aK (q'q). (8) 

^4/^4/4 

here, {t/t/t/"}^^ are the Racah coefficients in Eq. (jl]); {CCC}isf are the isoscalar factors 
in Eq. ([7]). All the transformation coefficients are tabulated in the appendix. 

A physical 5-quark state with quantum number a = iYIJ) is expressed as a channel 
coupling wave function 

^(^) = E^^(g 4 g). (9) 

k 

The channel coupling coefficients Ck is determined by the diagonalization of the 5-quark 

Hamiltonian as usual. As before, the calculation of Hamiltonian matrix elements in the 

cluster bases is replaced by that of matrix elements in the symmetry bases, where the 

fractional parentage technique can be used. The calculation of the matrix elements in 

I I 

symmetry bases is given in Ref. [5J, we will not repeat the procedure here. 



III. QUARK MODELS AND CALCULATION METHOD 

The models used in the calculations are naive quark model, Salamanca chiral quark model 
and quark derealization color screeningmodel. The Hamiltonian of these models have been 
given in the first paper of this series [5|. The only modification is the delocalized single 
particle orbital wavefunction, which is given below for the baryon-meson basis: 

^(r) = (0 L + ei fi )/A( ei ), 

Vv(r) = (0H + e 2 0L)/A(e 2 ), (10) 

N(e x ) = v/l + e? + 2e 1 ( ( p R \ ( p L ). (11) 

A(e 2 ) = v/l + ei + 2e 2 (0 fl |0 L ). 

where e\ denotes the quark in baryon orbit delocalized into meson orbit, while e 2 is from 
meson to baryon. They are determined by the dynamics of the 5-quark system. 
The calculation method is quite the same as in Ref. j^]. 

5 



IV. RESULTS AND DISCUSSIONS 



To take into consideration of the effects of various color structures, the hidden color 
channels are included in this calculation, which is different from our dibaryon calculation 
where only color singlet channels are included % [l2|. Like the systematic study of the JW 
di-quark configuration, all possible states within the u, d, s three-flavor world have been 
calculated in this baryon-meson configuration. Both single-channel and channel coupling 
calculations have been carried out in three quark models. To save space, only several general 
features and comparison between the JW di-quark and meson-baryon two configurations and 
among three models are given below. 

(1) Generally there exist effective attractions for almost all the states both in the extended 
QDCSM and the chiral quark model, while there are about ten channels which are pure 
repulsive in the naive quark model (NN interaction study already showed that the naive 
quark model is not realistic for hadron interaction, it can not provide the intermediate range 
NN attraction). In QDCSM, the attraction between decuplet baryons and vector mesons 
are usually large (> 200 MeV), the attraction between octet baryons and pseudo-scalar 
mesons are small (less than 10 MeV), and the attraction between octet (decuplet) baryons 
and vector (pseudo-scalar) mesons lie between. There are also several exceptions. This 
situation is similar to the baryon-baryon interactions where the attraction between decuplet 
baryons is large, between octet baryons is small and between decuplet and octet baryons 
is in between. In chiral quark model, the above general features are kept but with more 
exceptions. Fig.l gives several examples. It is clear that Ap,AK* and Yi*K* channels have 
strong attractions, whereas Utt, YiK and "EK show pure repulsive or very small attractions. 
More examples can be found in Fig. 3-6. 
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Fig.l Effective potentials vs cluster separation in QDCSM. 

(2) For chiral quark model and naive quark model, most states have very similar results 
(energy difference less than 30 MeV) for two configurations: qq-qq-q and qqq-qq. For these 
states, the separations between quark clusters corresponding to the minimum-energy are 
small (S and/or T < 0.7 fm). The differences between two configurations tend to disappear 
after channel-coupling, the remaining minor differences come from the different quark orbital 
wave functions used. There are also several states, the energy differences are large (greater 
than 60 MeV, even 100 MeV) for two configurations. For these states, the separation between 
quark clusters corresponding to minimum-energy in baryon-meson configuration is a little 
larger (color-singlet clusters can stay far away), while the ones in JW di-quark configuration 
are smaller (colorful clusters tend to stay close), that makes the energies different. However, 
the extended QDCSM gives different results. The di-quark configuration always gives lower 
energy than baryon-meson configuration. The main reason is that the color confinement 
is screened between quark pair in different clusters in QDCSM. The more the number of 
clusters, the lower the model energy. The screened color confinement assumed for quark 
pair in two color singlet clusters was shown to be reasonable by the previous studies 4, f|, 
however extends to colorful clusters may be questionable. One may expect there might be 
additional energy for separating color clusters. Further investigation is needed. 

(3) For the states in anti-decuplet, 27-plet and 35-plet the mass hierarchies obtained in 
meson-baryon configuration are quite similar to that of JW di-quark configuration. Here 
only the mass hierarchies of 8/ + 10/ states with J p = | + of two configurations in three 



quark models are given in Fig. 2. 
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Fig.2 mass hierarchy in the 8/ + 10/ with J 



in three quark models, compared with 



the one in diquark-diquark configuration and Ref. 
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(4) Figs. 3-6 show the effect of channel coupling (sc stands for single channel, cs for color 
singlet channel coupling and cc for full channel coupling). For most channels, the effect of 
channel coupling is small, especially for QDCSM. It is worth to note that the hidden-color 
states have minor effect on the lowest energy of a channel. Figs. 7-8 show that channel 
coupling effect for (YIJ) = (00|), where the channel coupling contributes a considerable 
attraction. However the contribution of hidden-color channels are still negligible, especially 
for QDCSM. We like to emphasize that this result is based on the assumption that the quark 
interaction used for hadron spectroscopy and color singlet hadron channels can be directly 
extend to hidden color channels and this assumption is questionable. In fact we don't have 
any information about the quark interaction between hidden color channels and colorless 
ones up to now. 
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Fig.3 Effective potentials for (YIJ) = (20±). 
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Fig.4 Effective potentials for (YIJ) = (2l|) in QDCSM. 
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Fig. 5 Same as fig.4 for chiral quark model. 
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Fig.6 Effective potentials (YIJ) = in QDCSM. 
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Fig.7 Effective potentials for (YIJ) = (00|) in QDCSM. 
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Fig.8 Effective potentials for (YIJ) = (00|) in chiral quark model. 

(5) Generally there exist attractions for baryon-meson system, however in most cases, 
the attraction is not enough to make the energy lower than the threshold, i.e., the sum of 
corresponding baryon and meson mass. Therefore there will be no narrow resonances except 
there is special reason to prevent the decay. Take + for example, the lowest energy state 
here is the KN state. The effective potential for the N — K channels in three quark models 
is given in Fig. 3. For the 1 = S-wave KN state, although the naive quark model gives 
repulsive potentials, there is a moderate attraction ~ 51 MeV at a separation of 0.6 fm 
in the QDCSM and a little weaker attraction ~ 48 MeV at a separation of 0.6 fm in the 
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chiral quark model. (The results are different from the previous calculation [141 ]. where no 



Q, 



attraction was obtained. The reason is the effect of a exchange was reduced in Ref. [14j by 
using larger mass and ideal mixing of o§ and erg). In QDCSM, the mass of B + is 1851 MeV, 
however the mass of kaon is poorly estimated to be 647 MeV in QDCSM. If we make a 
correction of the physical mass of kaon, the energy of the state can be reduced to 1700 MeV, 
which is still too high to match the experimental value 1540 MeV. In the chiral quark model, 
the case is quite similar to QDCSM. 

(6) Finally, there are several interesting states in QDCSM: (1) Their masses are lower 
than the lowest two- and three-body decay threshold and they might be better pentaquark 
candidates. However the results are sensitive to model details. For (YIJ) = (l^f ) channel, 
the lowest state has mass around 1822 MeV, which is not only lower than the two-body 
decay thresholds (Ap), but also lower than the three-body decay thresholds (Nup). And 
for (YIJ) = (00 1 ) channel, there is a state with mass around 2017 MeV, which is lower 
than the corresponding two-body decay threshold S*p and the three-body decay thresholds 
(A717?). (2) K~p state has been calculated dynamically in order to see if there is a resonance 
state A (1405). The mass of K~p is 1560 MeV. If we make a correction of the physical mass 
of kaon, the energy of the state can be reduced to 1408 MeV, which corresponds to the 
observed mass of A(1405). (3) In Ref. [151], several quasibound states have been proposed , 
such as 1=1 AK, I = | T,K, I — ^, § Ncj), and so on. In QDCSM, all these states have been 
calculated and strong attraction (> 100 MeV) were obtained. The dynamical calculation of 
these states are in progress. 



V. SUMMARY 



Multi-quark system is complicate, our knowledge about multi-quark system is limited 
and there is no well established multi-quark state experimentally up to now. Nevertheless 
its study is indispensable if one wants to understand the properties of the abundant color 
structures allowed in QCD. The systematic study could provide the general features of the 
multi-quark system. In this case, a powerful method is needed. This paper (and the previous 
onefsl) reports a group theoretic method for five-quark calculation. The method is applied 
to a systematic study of all possible five-quark states within u, d, s three-flavor world in 
baryon-meson configuration with three quark models. The powerful feature of group theory 
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method is shown by the large amount of spectroscopy data obtained easily in this approach 
(only very limited results have been shown in this report). 

In the pentaquark studies few years ago various constituent quark models were proposed. 
Our transformation between different physical bases to the same set of symmetry bases 
shows these models are a choice of different parts of the model space of the five-quark 
system. Our calculation shows that if one uses a large enough model space the differences 
of these different configuration choices tend to disappear. 

Besides the five-quark calculation, the baryon-meson scattering is another field where 
the technique developed in these two papers can be used. With the progress on the meson 
spectroscopy of non-relativistic quark model [if]], the unified description of hadron properties 
and hadron-hadron interaction becomes possible. There are a lot of work devoted to baryon- 
baryon scattering, baryon-meson scattering is less studied. Employing the group theoretic 
method presented here one can study the baryon-meson scattering within the framework of 
constituent quark models. 

There are experimental indication and theoretical motivation that the pure q 3 baryon and 
qq meson configuration should be modified to include higher Fock components. The effect 
of the five-quark component q 3 qq of baryon can be studied as well by the same technique 
developed here. 

This work is supported by the NSFC 10375030, 90503011, 10435080, 10775072. 
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Appendix 

Table Al. Transformation coefficients between physical bases and symmetry bases. The 
column labels are [z/ 4 ], [// 4 ], [/ 4 ], [oj 4 ], [J 4 ]. For first four labels, 1 stands for the symmetry 
[4]; 2, [31]; 3, [22]; 4, [211], and for the last one, 1 stands for isospin 2; 2, §; 3, 1; 4, \; 5, 
0. The row labels are B (baryon index) and M (meson index), where indices B and M are 
listed in Table A2. The transformation coefficients should be the square root of the entries, 
an d a negative sign mea ns to take the negative square root. 
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Table A2. Indices of baryons and mesons. The subscript 8 means color octet; the super- 



script | refers to the spin of baryon. A s is the flavor singlet A. 
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